Stein's method meets Malliavin calculus: 
a short survey with new estimates 

by Ivan NourdirE] and Giovanni Peccatlll 

Universite Paris VI and Universite Paris Quest 

Abstract: We provide an overview of some recent techniques involving the MalHavin calculus 
of variations and the so-called "Stein's method" for the Gaussian approximations of probability 
distributions. Special attention is devoted to establishing explicit connections with the classic 
method of moments: in particular, we use interpolation techniques in order to deduce some new 
estimates for the moments of random variables belonging to a fixed Wiener chaos. As an illus- 
tration, a class of central limit theorems associated with the quadratic variation of a fractional 
Brownian motion is studied in detail. 

Key words: Central limit theorems; Fractional Brownian motion; Isonormal Gaussian processes; 
Malliavin calculus; Multiple integrals; Stein's method. 

2000 Mathematics Subject Classification: 60F05; 60G15; 60H05; 60H07. 



Contents 

1 Introduction 

1.1 Stein's heuristic and method 

1.2 The role of Malliavin calculus . . . . 

1.3 Beyond the method of moments . . . 

1.4 An overview of the existing literature 



Preliminaries 

2.1 Isonormal Gaussian processes 

2.2 Chaos, hypercontractivity and products 

2.3 The language of Malliavin calculus . . 



One-dimensional approximations El 

3.1 Stein's lemma for normal approximations 11 

3.2 General bounds on the Kolmogorov distance 13 

3.3 Wiener chaos and the fourth moment condition 14 

3.4 Quadratic variation of the fractional Brownian motion, part one 17 

3.5 The method of (fourth) moments: explicit estimates via interpolation ... l20 



*Laboratoire de Probabilites et Modeles Aleatoires, Universite Pierre et Marie Curie, Boite courrier 
188, 4 Place Jussieu, 75252 Paris Cedex 5, France, ivan.nourdinOupmc.fr 

^Equipe ModaFX, Universite Paris Ouest - Nanterre la Defense, 200 Avenue de la Republique, 92000 
Nanterre, and LSTA, Universite Paris VI, France. Email: giovanni .peccatiOgmail . com 



1 



Multidimensional case u 

4.1 Main bounds 23 

4.2 Quadratic variation of fractional Brownian motion, continued l26 



1 Introduction 

This survey deals with the powerful interaction of two probabilistic techniques, namely 
the Stein's method for the normal approximation of probability distributions, and the 
Malliavin calculus of variations. We will first provide an intuitive discussion of the theory, 
as well as an overview of the literature developed so far. 

1.1 Stein's heuristic and method 

We start with an introduction to Stein's method based on moments computations. Let 
~ ^{0, 1) be a standard Gaussian random variable. It is well-known that the (integer) 
moments of N, noted fip := E{Np) for p ^ 1, are given by: /ij, = if p is odd, and 
f^p = {p ~ 1)" •= pV(2^''^(p/2)!) if p is even. A little inspection reveals that the sequence 
{/ip : p ^ 1} is indeed completely determined by the recurrence relation: 

/ii = 0, /i2 = 1, and fip = {p — 1) X /ip_2, for every p ^ 3. (1.1) 

Now (for p ^ 0) introduce the notation fp{x) = x^, so that it is immediate that the relation 
fll.ip can be restated as 

E[N X U-i{N)] = E[f;^,{N)], for every p^l. (1.2) 

By using a standard argument based on polynomial approximations, one can easily prove 
that relation (11.21) continues to hold if one replaces fp with a sufficiently smooth function 
/ (e.g. any function with a sub-polynomial derivative will do). Now observe that a 
random variable Z verifying E[Z fp_i{Z)] = E[fp_^{Z)] for every p ^ 1 is necessarily such 
that E{ZP) = Hp for every p ^ 1. Also, recall that the law of a o/K(0, 1) random vari- 
able is uniquely determined by its moments. By combining these facts with the previous 
discussion, one obtains the following characterization of the (standard) normal distribu- 
tion, which is universally known as "Stein's Lemma": a random variable Z has a ^(0, 1) 
distribution if and only if 

E[Zf{Z)-f'{Z)]=0, (1.3) 

for every smooth function f. Of course, one needs to better specify the notion of "smooth 
function" - a rigorous statement and a rigorous proof of Stein's Lemma are provided at 
Point 3 of Lemma [3. II below. 

A far-reaching idea developed by Stein (starting from the seminal paper [36|) is the fol- 
lowing: in view of Stein's Lemma and given a generic random variable Z, one can measure 
the distance between the laws of Z and A^ ~ -/^(0, 1), by assessing the distance from zero 
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of the quantity E[Zf{Z) — f'{Z)]^ for every / belonging to a "sufficiently large" class of 
smooth functions. Rather surprisingly, this somewhat heuristic approach to probabilistic 
approximations can be made rigorous by using ordinary differential equations. Indeed, one 
of the main findings of [36] and [37] is that bounds of the following type hold in great 
generality: 



where: (i) Z is a generic random variable, (ii) ~ ^(0,1), (iii) d{Z,N) indicates an 
appropriate distance between the laws of Z and (for instance, the Kolmogorov, or the 
total variation distance), (iv) JF is some appropriate class of smooth functions, and (v) C 
is a universal constant. The case where d is equal to the Kolmogorov distance, noted dxoi, 
is worked out in detail in the forthcoming Section 13. It we anticipate that, in this case, 
one can take C = 1, and JF equal to the collection of all bounded Lipschitz functions with 
Lipschitz constant less or equal to 1. 

Of course, the crucial issue in order to put Stein-type bounds into effective use, is 
how to assess quantities having the form of the RHS of (11.41) . In the last thirty years, 
an impressive panoply of approaches has been developed in this direction: the reader is 
referred to the two surveys by Chen and Shao [7] and Reinert [33] for a detailed discussion 
of these contributions. In this paper, we shall illustrate how one can effectively estimate 
a quantity such as the RHS of (11.41) . whenever the random variable Z can be represented 
as a regular functional of a generic and possibly infinite-dimensional Gaussian field. Here, 
the correct notion of regularity is related to Malliavin-type operators. 

1.2 The role of Malliavin calculus 

All the definitions concerning Gaussian analysis and Malliavin calculus used in the Intro- 
duction will be detailed in the subsequent Section [2l Let X = {X{h) : /i e i^} be an 
isonormal Gaussian process over some real separable Hilbert space Sj. Suppose Z is a 
centered functional of X, such that E{Z) = and Z is differentiable in the sense of Malli- 
avin calculus. According to the Stein-type bound ( 11.41) . in order to evaluate the distance 
between the law of Z and the law of a Gaussian random variable A^ ~ ^(0, 1), one must 
be able to assess the distance between the two quantities E[Zf{Z)] and E[f'{Z)]. The 
main idea developed in [TH], and later in the references [HI [20l [22l [23] , is that the needed 
estimate can be realized by using the following consequence of the integration by parts 
formula of Malliavin calculus: for every / sufficiently smooth (see Section 12.31 for a more 
precise statement), 



where D is the Malliavin derivative operator, is the pseudo-inverse of the Ornstein- 
Uhlenbeck generator, and (■, is the inner product of ^. It follows from (11.51) that, if the 
derivative /' is bounded, then the distance between E[Zf{Z)] and E[f'{Z)] is controlled by 




(1.4) 



E[Zf{Z)] = E[fiZ){DZ,-DL-'Z)^] 



(1.5) 
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the L^(r2)-norm of the random variable 1 — {DZ, —DL^^Z)^,. For instance, in the case of 
the Kolmogorov distance, one obtains that, for every centered and Malliavin differentiable 
random variable Z, 



We will see in Section [3731 that, in the particular case where Z = Iq{f) is a multiple Wiener- 
Ito integral of order g ^ 2 (that is, Z is an element of the gth Wiener chaos of X) with 
unit variance, relation HIM yields the neat estimate 



Note that E{Z^) — 3 is just the fourth cumulant of Z, and that the fourth cumulant of 
equals zero. We will also show that the combination of (11. 6p and Hl.lfi allows to recover 
(and refine) several characterizations of CLTs on a fixed Wiener chaos - as recently proved 
in [26] and [27]. 

1.3 Beyond the method of moments 

The estimate (11.71) . specially when combined with the findings of [23] and [3l] (see Section 
S]), can be seen as a drastic simplification of the so-called "method of moments and cumu- 
lants" (see Major [13] for a classic discussion of this method in the framework of Gaussian 
analysis). Indeed, such a relation implies that, if {Zn : n ^ 1} is a sequence of random 
variables with unit variance belonging to a fixed Wiener chaos, then, in order to prove 
that Zn converges in law to ~ ^(0, 1), it is sufficient to show that E{Z^) converges to 
E^N"^) = 3. Again by virtue of I11.7L one also has that the rate of convergence of E{Z^) to 
3 determines the "global" rate convergence in the Kolmogorov distance. In order to further 
characterize the connections between our techniques and moments computations, in Propo- 
sition [SHH we will deduce some new estimates, implying that (for Z with unit variance and 
belonging to a fixed Wiener chaos), for every integer k ^ 3 the quantity \E[Z^] — E[N^]\ 
is controlled (up to an explicit universal multiplicative constant) by the square root of 
\E[Z^] — E[N^]\. This result is obtained by means of an interpolation technique, recently 
used in [22] and originally introduced by Talagrand - see e.g. [38] . 

1.4 An overview of the existing literature 

The present survey is mostly based on the three references [I8], [22] and [23], dealing with 
upper bounds in the one-dimensional and multi-dimensional approximations of regular 
functionals of general Gaussian fields (strictly speaking, the papers [18] and [22] also contain 
results on non-normal approximations, related e.g. to the Gamma law). However, since 
the appearance of [H], several related works have been written, which we shall now shortly 
describe. 



dKoi{Z, N)i^E\l- {DZ, -DL^'Z)^ 



(1.6) 




(1.7) 
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- Our paper [19] is again based on Stein's method and Malliavin calculus, and deals 
with the problem of determining optimal rates of convergence. Some results bear 
connections with one-term Edgeworth expansions. 

- The paper [3], by Breton and Nourdin, completes the analysis initiated in [HI Section 
4], concerning the obtention of Berry-Esseen bounds associated with the so-called 
Breuer-Major limit theorems (see [5]). The case of non-Gaussian limit laws (of the 
Rosenblatt type) is also analyzed. 

- In [20], by Nourdin, Peccati and Reinert, one can find an application of Stein's 
method and Malliavin calculus to the derivation of second order Poincare inequalities 
on Wiener space. This also refines the CLTs on Wiener chaos proved in [26] and [27] . 

- One should also mention our paper [16], where we deduce a characterization of non- 
central limit theorems (associated with Gamma laws) on Wiener chaos. The main 
findings of [H] are refined in [T8] and [22], again by means of Stein's method. 

- The work [24], by Nourdin and Viens, contains an application of (11.51) to the estimate 
of densities and tail probabilities associated with functionals of Gaussian processes, 
like for instance quadratic functionals or suprema of continuous-time Gaussian pro- 
cesses on a finite interval. 

- The findings of [24] have been further refined by Viens in [40], where one can also 
find some applications to polymer fiuctuation exponents. 

- The paper [4], by Breton, Nourdin and Peccati, contains some statistical applica- 
tions of the results of [24], to the construction of confidence intervals for the Hurst 
parameter of a fractional Brownian motion. 

- Reference [2], by Bercu, Nourdin and Taqqu, contains some applications of the results 
of [H] to almost sure CLTs. 

- In [21], by Nourdin, Peccati and Reinert, one can find an extension of the ideas 
introduced in [18] to the framework of functionals of Rademacher sequences. To this 
end, one must use a discrete version of Malliavin calculus (see Privault [32]). 

- Reference [29], by Peccati, Sole, Taqqu and Utzet, concerns a combination of Stein's 
method with a version of Malliavin calculus on the Poisson space (as developed by 
Nualart and Vives in [28]). 

- Reference [22], by Nourdin, Peccati and Reinert, contains an application of Stein's 
method, Malliavin calculus and the "Lindeberg invariance principle", to the study 
of universality results for sequences of homogenous sums associated with general 
collections of independent random variables. 
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2 Preliminaries 



We shall now present the basic elements of Gaussian analysis and Malliavin calculus that 
are used in this paper. The reader is referred to the monograph by Nualart [25] for any 
unexplained definition or result. 

2.1 Isonormal Gaussian processes 

Let ^ be a real separable Hilbert space. For any g ^ 1, we denote by S)'^'^ the gth 
tensor product of S), and by the associated gth symmetric tensor product; plainly, 

We write X = {X{h), h E S)} to indicate an isonormal Gaussian process over S). This 
means that X is a centered Gaussian family, defined on some probability space JF, P), 
and such that E [X{g)X{h)] = {g,h)f, for every g,h E Sj. Without loss of generality, we 
also assume that is generated by X. 

The concept of an isonormal Gaussian process dates back to Dudley's paper [10]. As 
shown in the forthcoming five examples, this general notion may be used to encode the 
structure of many remarkable Gaussian families. 

Example 2.1 (Euclidean spaces) Fix an integer d ^ 1, set = W'- and let (ei, ...,6^) 
be an orthonormal basis of R"' (with respect to the usual Euclidean inner product). Let 
{Zi,...,Zct) be a Gaussian vector whose components are i.i.d. A^(0, 1). For every h = 
Yl'j=i ^i^j (where the Cj are real and uniquely defined), set X [h) = Yl'j=i ^j^j define 
X = [X (h) : h eW^Y Then, X is an isonormal Gaussian process over M'^ endowed with 
its canonical inner product. 

Example 2.2 (Gaussian measures) Let {A,A,i') be a measure space, where z/ is posi- 
tive, (T-finite and non-atomic. Recall that a (real) Gaussian random measure over {A, A), 
with control is a centered Gaussian family of the type 

a = {G{B) : B e A, v{B) < oo} , 

satisfying the relation: for every B,C G ^ of finite z^-measure, E[G{B)G{C)] = u{B fl 
C). Now consider the Hilbert space = L"^ {A,A,h'), with inner product {h,h')sj = 
f^h{a)h'{a)h'{da). For every h E Sj, define X (h) = f^h{a)G{da) to be the Wiener-Ito 
integral of h with respect to G. Then, X = {X (h) : h E (Z, Z, u)} defines a centered 
Gaussian family with covariance given by E[X{h)X{h')] = {h,h')fj, thus yielding that X 
is an isonormal Gaussian process over {A, A, p). For instance, by setting A = [0, +oo) 
and z/ equal to the Lebesgue measure, one obtains that the process Wt = G{[0,t)), t ^ 0, 
is a standard Brownian motion started from zero (of course, in order to meet the usual 
definition of a Brownian motion, one has also to select a continuous version of W), and X 
coincides with the L^(f2)-closed linear Gaussian space generated by W. 
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Example 2.3 (Isonormal spaces derived from covariances) Let Y = {Yt : t ^ 0} be 

a real-valued centered Gaussian process indexed by the positive axis, and set R{s,t) = 
E \YsYt\ to be the covariance function of Y . One can embed Y into some isonormal Gaus- 
sian process as follows: (i) define S" as the collection of all finite linear combinations of 
indicator functions of the type l[o,t], t ^ 0; (ii) define = SjRtohe the Hilbert space given 
by the closure of ^ with respect to the inner product 

(/, h)R := ^aiCjR{si,tj) , 

where / = ^jflilp.si] and h = ^jCjl[o^tj] are two generic elements of tf; (iii) for h = 
^jCjl[o^tj] e ^, set X (h) = Y^jCjYtr^ (iv) for h G ^r, set X {h) to be the (P) limit 
of any sequence of the type X{hn), where {/i„} C converges to h in S)r. Note that 
such a sequence {/i„} necessarily exists and may not be unique (however, the definition of 
X {h) does not depend on the choice of the sequence {hn}). Then, by construction, the 
Gaussian space {X (h) : h E S^r} is an isonormal Gaussian process over S^r. See Janson 
[T2l Ch. 1] or Nualart [25], as well as the forthcoming Section [TH for more details on this 
construction. 

Example 2.4 (Even functions and symmetric measures) Other classic examples of 
isonormal Gaussian processes (see e.g., [6l[IIl[l3]) are given by objects of the type 

x^ = {Xf, W-.ije s^eA , 

where /3 is a real non-atomic symmetric measure on (— 7r,7r] (that is, j3 {dx) = [3{—dx)), 
and 

9)E,p = Ll{{-'K,Tr],df3) (2.8) 

stands for the collection of all real linear combinations of complex-valued even functions 
that are square-integrable with respect to (3 (recall that a function ip is even if ip {x) = 
ip {—x)). The class S^e,i3 is a real Hilbert space, endowed with the inner product 

^i(x)7/'2(-x)/3(rfx) gM. (2.9) 

-vr 

This type of construction is used in the spectral theory of time series. 

Example 2.5 (Gaussian Free Fields) Let ^ 2 and let D be a domain in R'^. Denote 
by Hs{D) the space of real- valued continuous and continuously differentiable functions on 
M"^ that are supported on a compact subset of D (note that this implies that the first 
derivatives of the elements of Hs{D) are square-integrable with respect to the Lebesgue 
measure). Write H{D) in order to indicate real Hilbert space obtained as the closure of 
Hs{D) with respect to the inner product (/,(?) = /j^^ V/(x) ■ Vg{x)dx^ where V is the 
gradient. An isonormal Gaussian process of the type X = {X{h) : h G H{D)} is called 
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a Gaussian Free Field (GFF). The reader is referred to the survey by Sheffield [35] for a 
discussion of the emergence of GFFs in several areas of modern probability. See e.g. Rider 
and Virag [34] for a connection with the "circular law" for Gaussian non-Hermitian random 
matrices. 

Remark 2.6 An isonormal Gaussian process is simply an isomorphism between a centered 
L^(r2)-closed linear Gaussian space and a real separable Hilbert space S}. Now, fix a generic 
centered L^(f2)-closed linear Gaussian space, say Q. Since Q is itself a real separable Hilbert 
space (with respect to the usual L^(fi) inner product) it follows that Q can always be 
(trivially) represented as an isonormal Gaussian process, by setting = Q. Plainly, the 
subtlety in the use of isonormal Gaussian processes is that one has to select an isomorphism 
that is well-adapted to the specific problem one wants to tackle. 



2.2 Chaos, hypercontractivity and products 

We now fix a generic isonormal Gaussian process X = {X{h),h G defined on some 
space (f2, JF, P) such that cr(X) = JF. 

Wiener chaos. For every g ^ 1, we write Hg in order to indicate the gth Wiener chaos of 
X. We recall that Tig is the closed linear subspace of L^(f2, JF, P) generated by the random 
variables of the type Hq{X{h)), where h E is such that \\h\\^ = 1, and Hg stands for the 
gth Hermite polynomial, defined as 

Hg{x) = {-lye^^e-"^, xER, q ^ 1. (2.10) 

We also use the convention Hq = M. For any g ^ 1, the mapping 

Ig{h^'^)=q\Hg{X{h)) (2.11) 

can be extended to a linear isometry between the symmetric tensor product i^®'^ equipped 
with the modified norm v^lMlf,®'? ^ind the gth Wiener chaos Hg. For g = 0, we write 
/o(c) =c,ceR. 

Remark 2.7 When = L"^ (A, A, z/), the symmetric tensor product i^®"^ can be identified 
with the Hilbert space {A'^^A'^, u'^), which is defined as the collection of all symmetric 
functions on A'^ that are square-integrable with respect to u'^. In this case, it is well-known 
that the random variable Ig{h), h G ^®'^, coincides with the (multiple) Wiener-Ito integral, 
of order q, of h with respect to the Gaussian measure B i— X{1b), where B E A has finite 
i^-measure. See [25l Chapter 1] for more details on this point. 

Hypercontractivity. Random variables living in a fixed Wiener chaos are hypercontrac- 
tive. More precisely, assume that Z belongs to the gth Wiener chaos Hg (g ^ 1). Then, 
Z has a finite variance by construction and, for all p G [2,+oo), one has the following 
estimate (see [121 Th. 5.10] for a proof): 

E{\Z\P) ^{p- 1)P1/^E{ZY^^- (2.12) 
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In particular, if E{Z^) = 1, one has that E{\Z\p) ^ (p - l)^"'/^. por future use, we also 
observe that, for every g ^ 1, the mapping p t-^ {p— 1)P'J/2 jg strictly increasing on [2, +oo). 

Chaotic decompositions. It is well-known {Wiener chaos decomposition) that the space 
JF, P) can be decomposed into the infinite orthogonal sum of the spaces Hg. It follows 
that any square-integrable random variable Z G L'^{Q, J^, P) admits the following chaotic 
expansion 

oo 

Z = Y.Wg), (2.13) 

where /o = E{Z), and the kernels fq G ^l®'', q ^ I, are uniquely determined. For every 
g ^ 0, we also denote by Jq the orthogonal projection operator on Tiq. In particular, if 
Z E L'^{n,J^,P) is as in (l2:T3D . then Jq{Z) = Iq{fq) for every g ^ 0. 

Contractions. Let {e^. A; ^ 1} be a complete orthonormal system in i^. Given / G 
and g G for every r = 0, . . . , p A g, the contraction of / and g of order r is the element 
of ^®(p+<?-2r) defined by 

oo 

f^r9= 5Z (/,eii ® . . . ® e,Ji=,®r ® (5(,ei, ® . . . ® e,Ji3«r. (2.14) 

jl,...,ir = l 

Notice that / ®r g is not necessarily symmetric: we denote its symmetrization by f^rd G 
^0(p+g-2r) Moreover, f ®o9 = f®9 equals the tensor product of / and g while, for p = g, 
one has that f ®q g = {f,g)f,^i- In the particular case where = L^{A,A,i'), one has 
that ^®'^ = Ll{A'^ , A'^ , u'^) (see Remark \2.1\ \ and the contraction in (12.141) can be written 
in integral form as 



(/ ®r fi')(^l5 • • • 5 ^p+q-2r) — / /(^l, • • • 5 ^p-r, ^i, . . . , S^) 

X 5'(^p-r+l; • • • ; tp+q-2r, ^l, . . . , Sr)dh'{Si) . . . dv{Sr). 



Multiplication. The following multiplication formula is well-known: if / G ii®^ and 
g G then 

W)k{9) = Ilj\r) (^J W2.(/®.^7). (2.15) 

Note that (I2.15P gives an immediate proof of the fact that multiple Wiener-Ito integrals 
have finite moments of every order. 

2.3 The language of Malliavin calculus 

We now introduce some basic elements of the Malliavin calculus with respect to the isonor- 
mal Gaussian process X. 
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Malliavin derivatives. Let S be the set of all cylindrical random variables of the type 



Z=r7(X(0i),...,X(0O), (2.16) 

where ^ 1, (7 : R" ^ M is an infinitely differentiable function with compact support and 
(f)i G io. The Malliavin derivative of Z with respect to X is the element of L^(fi, Sj) defined 
as 

1=1 * 

By iteration, one can define the mth derivative D"^Z, which is an element of 

for every m ^ 2. For m ^ 1 and p ^ 1, D™'^ denotes the closure of S with respect to the 

norm || ■ \\m,p, defined by the relation 



\Zrm,p = E[\Zn+J2E{\\D^Z\\l,^). 



1=1 



The chain rule. The Malliavin derivative D verifies the following chain rule. If ip : 
M is continuously differentiable with bounded partial derivatives and if Z = {Zi, . . . , Za) 
is a vector of elements of 3^'^, then (p{Z) G D^'^ and 

Dv(^) = i2^(^)DZ,. (2.17) 

i=l * 

A careful application e.g. of the multiplication formula f l2.15p shows that f l2.17p continues 
to hold whenever the function (/? is a polynomial in d variables. Note also that a random 
variable Z as in (I2.13P is in D^'^ if and only if 9||</q(2^)||^2(Q) < 00 and, in this 

case, E(||DZ|||) = X]^i ?ll'A(^)llL2(n)- If ^ = L'^{A,A,u) (with u non-atomic), then 
the derivative of a random variable Z as in (12.130 can be identified with the element of 
L'^{A X f2) given by 

00 

D,.Z = ^gVi (/,(-, x)), xeA. (2.18) 

9=1 

The divergence operator. We denote by S the adjoint of the operator D, also called the 
divergence operator. A random element u G L'^{Q,S^) belongs to the domain of S, noted 
Dom6, if and only if it verifies \E{DZ,u)sj\ ^ Cu \\Z\\L2(^fi^ for any Z G D^'^, where Cu is a 
constant depending only on u. If m G DomS, then the random variable S{u) is defined by 
the duality relationship 

E{Z6{u))=E{{DZ,u)s,), (2.19) 

which holds for every Z G D^'^. 
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Ornstein-Uhlenbeck operators. The operator L, known as the generator of the Ornstein- 

jq=0 



Uhlenbeck semigroup, is defined as L = J2T=o "l-^q- The domain of L is 



9=1 



There is an important relation between the operators D, 6 and L (see e.g. [25l Proposition 
1.4.3]): a random variable Z belongs to ©^'^ if and only if Z E Dom (5-0) (i.e. Z G ©^'^ 
and DZ G Dom5) and, in this case, 

6DZ = -LZ. (2.20) 

For any Z G Iv^(fi), we define L~^Z = J2'^=i~qJq{^)- The operator is called the 
pseudo-inverse of L. For any Z G iv^(f^), we have that L~^Z G DomL, and 

LL~^Z = Z-E(Z). (2.21) 

An important string of identities. Finally, let us mention a chain of identities playing 
a crucial role in the sequel. Let / : M — M be a function with bounded derivative, and 
let F, Z E Assume moreover that E{Z) = 0. By using successively fl2.2ip . fl2.20p and 
fl2.17p . one deduces that 

E{ZfiF)) = E{LL-'ZxfiF))=E{5Di-L-'Z)xf{F)) 
= E{{Df{F),-DL~'Z)^) 

= E{f{F){DF,-DL~'Z)^). (2.22) 

We will shortly see that the fact E{Zf{F)) = E{f{F){DF,-DL~^Z)^) constitutes a 
fundamental element in the connection between Malliavin calculus and Stein's method. 



3 One-dimensional approximations 
3.1 Stein's lemma for normal approximations 

Originally introduced in the path-breaking paper [36], and then further developed in the 
monograph [37], Stein's method can be roughly described as a collection of probabilis- 
tic techniques, allowing to characterize the approximation of probability distributions by 
means of differential operators. As already pointed out in the Introduction, the two surveys 
[7] and [33] provide a valuable introduction to this very active area of modern probability. 
In this section, we are mainly interested in the use of Stein's method for the normal approx- 
imation of the laws of real- valued random variables, where the approximation is performed 
with respect to the Kolmogorov distance. We recall that the Kolmogorov distance between 
the laws of two real- valued random variables Y and Z is defined by 

dKoi{Y, Z) = sup \P{Y <: z) - P{Z ^ z)\. 
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The reader is referred to [18] for several extensions of the results discussed in this 
survey to other distances between probability measures, such as e.g. the total variation 
distance, or the Wasserstein distance. The following statement, containing all the elements 
of Stein's method that are needed for our discussion, can be traced back to Stein's original 
contribution l36l. 



Lemma 3.1 Let N ~ ^(0, 1) be a standard Gaussian random variable. 
1. Fix z eM., and define : M ^ M as 

fz{x) = 6^^ f (l(_oo,.](«) - P{N ^ z))e-^da, x G M. (3.23) 



Then, fz is continuous onM, bounded by v27r/4, dijferentiable onM\{z}, and verifies 
moreover 

/:(x)-x/,(x) = l(_oo,.](x)-P(iV^2) foraUxER\{z}. (3.24) 
One has also that fz is Lipschitz, with Lipschitz constant less or equal to 1. 

2. Let Z be a generic random variable. Then, 

dKoi{Z,N) i^snv\E[Zf{Z) - f{Z)]\, (3.25) 
/ 

where the supremum is taken over the class of all Lipschitz functions that are bounded 
by -\/27r/4 and whose Lipschitz constant is less or equal to 1. 

3. Let Z be a generic random variable. Then, Z ~ ^(0, 1) if and only if E[Zf{Z) — 
f'{Z)] = for every continuous and piecewise dijferentiable function f verifying the 
relation E\f'{N)\ < oo. 

Proof. {Point 1) We shall only prove that fz is Lipschitz and we will evaluate its constant 
(the proof of the remaining properties is left to the reader). We have, for x ^ 0, x 7^ z: 



l(-oo,.](x) - P(Ar ^ z) + xe^ / (l(-oo,.](a) - P(iV ^ z))e-'^ da 

J —00 

2 2 
l(-oo,.](x) - P{N ^z)- xe^ / (l(-oo,.](a) - P(iV ^ z))e-^ da 



/ 2 2 

^ ||l(-oo,.i(-) - ^(^ ^ 2)|| (l+xe'^ e-^da 



X 



2 /-+°° 2 



^ 1 + 62 / ae 2 da = 2. 

J X 

Observe that identity (*) holds since 

I r+°° 2 
= E(l(_oo,.](iV) - PiN ^z))=-= (l(-oo,.](a) - P{N ^ z))e"-da. 

yj l-K J -00 
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For X ^ 0, X ^ z, we can write 



\fM\ 



X 



P{N ^ z)+xe' 



(1 



,1 a 



P{N ^ z))e''^da 



^ 1 



lale ^ da = 2. 



Hence, we have shown that fz is Lipschitz with Lipschitz constant bounded by 2. For the 
announced refinement (that is, the constant is bounded by 1), we refer the reader to Chen 
and Shao [7| Lemma 2.2]. 

{Point 2) Take expectations on both sides of (13.241) with respect to the law of Z. Then, 
take the supremum over all z G M, and exploit the properties of fz proved at Point 1. 

{Point 3) If Z ^(0, 1), a simple application of the Fubini theorem (or, equivalently, an 
integration by parts) yields that E[Zf{Z)] = E[f'{Z)] for every smooth /. Now suppose 
that E[Zf{Z) — f'{Z)] = for every function / as in the statement, so that this equality 
holds in particular for f = fz and for every 2; G M. By integrating both sides of fl3.24p 
with respect to the law of Z, this yields that P{Z ^ z) = P{N ^ z) for every 2; G M, and 
therefore that Z and N have the same law. ^ 



Remark 3.2 Formulae (13.241) and (I3.25P are known, respectively, as Stein's equation and 
Stein's bound. As already evoked in the Introduction, Point 3 in the statement of Lemma 
13. H is customarily referred to as Stein's lemma. 

3.2 General bounds on the Kolmogorov distance 

We now face the problem of establishing a bound on the normal approximation of a centered 
and Malliavin-differentiable random variable. The next statement contains one of the main 
findings of [T8] . 

Theorem 3.3 (See [H]) Let Z G D^'^ be such that E{Z) = and Var(Z) = 1. Then, 
forN ^ ^(0,1), 

dKoi{Z, N) ^ y/Var((DZ, -DL-^Z)^). (3.26) 

Proof. In view of (13.251) . it is enough to prove that, for every Lipschitz function / with 
Lipschitz constant less or equal to 1, one has that the quantity \E[Zf{Z) — f'{Z)]\ is less 
or equal to the RHS of (13.261) . Start by considering a function / : M ^ M which is and 
such that II /'I loo ^ 1. Relation (12:221) yields 

E{Zf{Z)) = E{f{Z){DZ,-DL-^Z)^), 
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so that 



\E{f{Z))-E{Zf{Z))\ = \E{f{Z){l-{DZ,-DL-'Z)s,))\ ^ E\1- {DZ, -DL-'Z)s,\. 

By a standard approximation argument (e.g. by using a convolution with an approxi- 
mation of the identity), one sees that the inequality — ^ E^l — 
{DZ, —DL~^Z)sj\ continues to hold when / is Lipschitz with constant less or equal to 1. 
Hence, by combining the previous estimates with (13.251) . we infer that 



dKoiiZ,N)^E\l-{DZ,-DL-^Z)^\ ^ ^JE{1- {DZ,-DL-^Z),,) . 

Finally, the desired conclusion follows by observing that, if one chooses f{z) = z in (12.221) . 
then one obtains 



E{{DZ, -DL~'Z),,) = EiZ') = l, 



(3.27) 



so that E 



(l - {DZ, -DL-^Z)^f = Var((DZ, -DL-^Z)sj). 



□ 



Remark 3.4 By using the standard properties of conditional expectations, one sees that 
(I3.26P also implies the "finer" bound 

dKoi{Z,N) ^ y/Var((7(Z)), (3.28) 

where g{Z) = E[{DZ, —DL^^Z)^\Z]. In general, it is quite difficult to obtain an explicit 
expression of the function g. However, if some crude estimates on g are available, then one 
can obtain explicit upper and lower bounds for the densities and the tail probabilities of 
the random variable Z. The reader is referred to Nourdin and Viens [24] and Viens [40] for 
several results in this direction, and to Breton et al. [1] for some statistical applications of 
these ideas. 



3.3 Wiener chaos and the fourth moment condition 

In this section, we will apply Theorem 13.31 to chaotic random variables, that is, random 
variables having the special form of multiple Wiener-Ito integrals of some fixed order g ^ 2. 
As announced in the Introduction, this allows to recover and refine some recent charac- 
terizations of CLTs on Wiener chaos (see [26l ET]). We begin with a technical lemma. 



Lemma 3.5 Fix an integer q ^ I, and let Z = Ig{f) (with f G S^®'^) be such that Var(Z) = 
E{Z'^) = 1. The following three identities are in order: 

^\\DZ\\l-l=qJ2ir-ly■('^Z^^ h,^2rU®rf), (3.29) 
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V8.T(^^\\DZ\\l^=''£j^r\'(^^^ (3.30) 



and 



E{Z') - 3 = -J^rrl' Q (2g - 2r)!||/® (3.31) 



r=l 

In particular, 



Var (^ll^^lll) ^ ^(^(^') - 3)- (3.32) 

Proof. Without loss of generality, we can assume that S) is equal to L'^{A,A,i'), where 
(^4, ^) is a measurable space and z/ a c-finite measure without atoms. For any a G A, we 
have DaZ = qlq-i{f{-,a)) so that 

-\\DZ\\l = Q I Vi(/(-,a))V(rfa) 

= g / Vrir /2,-2-2r(/(-,a)®,/(-,a))Kc?a) by (JSIS]) 

JAr=^ \ y 

= h,-2-2r(^jj{;a)®rf{;aHda)^ 



r=0 

9-1 / 1 \ 2 



r=0 ^ 

r=l ^ ^ 



g-1 

|2 



r=l ^ 



Since ^(Z^) = the proof of (13:291) is finished. The identity (ICTD follows from 

fl3.29p and the orthogonality properties of multiple stochastic integrals. Using (in order) 
formula 1^2^ and the relation D{Z'') = SZ^DZ, we infer that 

EiZ") = ^E{SDZ X Z^) = ^E{{DZ,D{Z'^))^) = ^E{Z^DZ\\l). (3.33) 

Moreover, the multiplication formula (12.151) yields 

= Wf = Y.s\ {f\ 'h,-2s{f ®s /). (3.34) 

By combining this last identity with (13.291) and ( 13.331) . we obtain (13.3ip and finally (13.321) . 
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□ 



As a consequence of Lemma 13.51 we deduce the following bound on the Kolmogorov 
distance - first proved in [22] . 

Theorem 3.6 (See |22| ) Let Z belong to the qth chaos Hq of X, for some q ^ 2. Suppose 
moreover that Var(Z) = E{Z'^) = 1. Then 

dKoiiZ, N) ^ ^^-^{E{Z^) (3.35) 

Proof. Since L'^Z = —\Z^ we have {DZ, —DL^^Z) fj = So, we only need to 

apply Theorem 1X3] and formula (13.321) . |--| 

The estimate (13.351) allows to deduce the following characterization of CLTs on Wiener 
chaos. Note that the equivalence of Point (i) and Point (ii) in the next statement was first 
proved by Nualart and Peccati in [27] (by completely different techniques based on stochas- 
tic time-changes), whereas the equivalence of Point (iii) was first obtained by Nualart and 
Ortiz-Latorre in [26] (by means of Malliavin calculus, but not of Stein's method). 

Theorem 3.7 (see |26L I27| ) Let (Z„) be a sequence of random variables belonging to the 
qth chaos Hg of X, for some fixed q ^ 2. Assume that Var(Z„) = E{Z^) = 1 for all n. 
Then, as n ^ oo, the following three assertions are equivalent: 

(i) Zn^Nr^ ^(0, 1); 

(ii) E{Zt) - E{N') = 3; 

(iii) Var(i||Z}Z„||2) ^0. 

Proof. For every n, write Zn = Iqifn) with /„ G uniquely determined. The implication 
(iii) — (i) is a direct application of Theorem l3.6l and of the fact that the topology of the Kol- 
mogorov distance is stronger than the topology of the convergence in law. The implication 
(i) (ii) comes from a bounded convergence argument (observe that sup^^j^ E{Z^) < oo 
by the hypercontractivity relation f l2.12p ). Finally, let us prove the implication (ii) (iii). 
Suppose that (ii) is in order. Then, by virtue of fl3.3ip . we have that ||/n®r/n|U®29-2r tends 
to zero, as n ^ oo, for all (fixed) r G — 1}. Hence, (13.301) allows to conclude that 

(iii) is in order. The proof of Theorem 13.71 is thus complete. i-. 



Remark 3.8 Theorem 13.71 has been applied to a variety of situations: see e.g. (but the 
list is by no means exhaustive) Barndorff-Nielsen et al. [I], Corcuera et al. [8], Marinucci 
and Peccati [H], Neuenkirch and Nourdin [I5], Nourdin and Peccati [I7| and Tudor and 
Viens [39], and the references therein. See Peccati and Taqqu [30] for several combinatorial 
interpretations of these results. 
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By combining Theorem 13.61 and Theorem 13.71 we obtain the following result. 

Corollary 3.9 Let the assumptions of Corollary \3. 71 prevail. As n ^ oo, the following 
assertions are equivalent: 

(a) Z„^iV~^(0,l); 

(b) dKoi{Zn,N)^0. 

Proof. Of course, only the implication (a) (b) has to be proved. Assume that (a) is in 
order. By Corollary 13.71 we have that Var — > 0. Using Theorem 13.61 we get 
that (b) holds, and the proof is done. |-| 



3.4 Quadratic variation of the fractional Brownian motion, part 
one 

In this section, we use Theorem l3.3l in order to derive an explicit bound for the second-order 
approximation of the quadratic variation of a fractional Brownian motion. 

Let B = {Bt : t ^ 0} be a fractional Brownian motion with Hurst index H G (0, 1). 
This means that i? is a centered Gaussian process, started from zero and with covariance 
function E^BgBt) = R{s,t) given by 

R{s, t) = i (t^^ + - |t - sp^) , s, t ^ 0. 

The fractional Brownian motion of index H is the only centered Gaussian processes nor- 
malized in such a way that Var(_Bi) = 1, and such that B is selfsimilar with index H and 
has stationary increments. If H = 1/2 then R{s,t) = min(s,t) and B is simply a standard 
Brownian motion. If if 7^ 1/2, then B is neither a (semi) martingale nor a Markov process 
(see e.g. [25] for more details). 

As already explained in the Introduction (see Example 12. 3p . for any choice of the 
Hurst parameter H ^ (0, 1) the Gaussian space generated by B can be identified with an 
isonormal Gaussian process X = {X{h) : h G S)}, where the real and separable Hilbert 
space Sj is defined as follows: (i) denote by the set of all M-valued step functions on 
[0, 00), (ii) define as the Hilbert space obtained by closing with respect to the scalar 
product 

(l[0,t], l[0,s])^ = R(t, S). 

In particular, with such a notation, one has that Bt = X(l[o,t]). 
Set 

1 'n?^ 

Zn = ^ [(^fc+l — — 1] =' ^ [{B(^k+\)/n — Bk/nf " 
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where cr„ > is chosen so that E{Zl) = 1. It is well-known (see e.g. [5]) that, for every 
H ^ 3/4 and for n — oo, one has that Z„ converges in law to ~ ^(0, 1). The following 
result uses Stein's method in order to obtain an explicit bound for the Kolmogorov distance 
between Z„ and A^. It was first proved in [H] (for the case H < 3/4) and [3] (for H = 3/4). 

Theorem 3.10 Let N ~ ^(0, 1) and assume that H ^ 3/4. Then, there exists a constant 
c/f > (depending only on H) such that, for every n ^ 1, 



dKol{Zn,N) ^ChX < 



n^^-i «/i/e[i,f) . (3.36) 

\/\og n ^ 4 



Remark 3.11 1- By inspection of the forthcoming proof of Theorem l3.10l one sees that 
lim^^oo ^ = 2 E.ez (r) HHe (0,3/4), with p given by I^M) , and lim„^oo ^ = 
9/16 if iJ = 3/4. 



2. When H > 3/4, the sequence (Z„) does not converge in law to ^{0, 1). Actually, 
Zn Zoo ~ "Hermite random variable" and, using a result by Davydov and Mar- 

n— >oo 

tynova [9], one can also associate a bound to this convergence. See [3j for details on 
this result. 

3. More generally, and using the analogous computations, one can associate bounds 
with the convergence of sequence 

^ n— 1 ^ n—1 

^n^ = -(^ X] ^Q^^k+l - Bk) = ^ Hq{n^{B(^k+l)/n " ^fc/n)) 

'^n k=o k=0 

towards A^ ~ ^(0,1), where Hg (g ^ 3) denotes the gth Hermite polynomial (as 
defined in (12.101) ). and an^ is some appropriate normalizing constant. In this case, 
the critical value is H = 1 — l/(2g) instead of H = 3/4. See [18] for details. 

In order to show Theorem 13. 101 we will need the following ancillary result, whose proof 
is obvious and left to the reader. 

Lemma 3.12 1. For r G Z, let 

p(r) = \{\r + Ip^ + |r - 1^^ - 2|r|2^). (3.37) 

If H ^ \, one has p{r) ~ H(2H — l)\r\'^^~^ as \r\ ^ oo. If H = \ and \r\ ^ 1, one 
has p{r) = 0. Consequently, X]rezP^(^) < oo i/ and only if H < 3/4. 

2. For all a > — 1, we have Ylr=i ~ n"'^^/ (a + 1) as n ^ oo. 
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We are now ready to prove the main result of this section. 
Proof of TheoremWM Since || |||, = E{{Bk+i - Btf) = 1, we have, by ^JTi), 



SO that Zn = hifn) with fn = -^ X]fc=o -'-^^fc+i] ^ Let us compute the exact value of 
(T„. Obse rve that (l[fc,fc+i], = E[{Bk+i - Bk){Bi+i - Bi)) = p{k - I) with p given 

by (13371) . Hence 



E 



'n-\ 



E 



,fc=0 

'n-\ \ 2 

.fc=0 



^ [(fifc+i - 5^)2 - 1] 

fc,/=0 
n-l 

2 5^ l[M+i])| = 2 5^ /(A; - /). 



[fc,fc+i]^ 



n-l 



That is, 



k,l=0 



n-l 



k,l=0 



71—1 n~l—l 



k,l=0 1=0 r=~l 

Assume that H < 3/4. Then, we have 

\r\ + 1~ 



2\n'£p\r)-Y,{\r\^l)p^ 

|r|<n |r|<n 



n 



L{|r-|<n}- 



Since X]rezP^(^) ^ obtain, by bounded Lebesgue convergence: 



lim ^ = 2 Vp2^r). 



(3.38) 



Assume that = 3/4. We have p^(r) ~ as |r| ^ oo. Therefore, as n — *■ cxo 



E2/ \ 9"- 1 On log 72 



|r|<n, 



0<|r|<7i 



32 



and 



|r'|<ri |r|<n 
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We deduce that 



lim 



9 



racoon log n 16 
Now, we have, see (13.301) for the first equahty, 

1 



Var(-pZ„||y 



2 1 /n ®1 fn\\%S)2 


1 


n-1 
k,l=0 


1 


n-1 
j,fc,/=0 


1 


n-1 
j,fc,Z=0 


^ n— 1 


2< 


/ n-1 
\s=-n+l 



1 

2^ 



n-1 

El 

fcj=0 



02 

[k,k+l] 



. 1 02 
■•l ^[1,1+1] 



n-1 



r=— n+1 



n-1 



r=— n+1 



(3.39) 



If ^ 1/2 then X]sezlP('^)l ^ J2r&zP^(^) < cx) so that, in view of (|3.38l) . 



Var(|||DZ„||2) = O(n-i). If 1/2 < iJ < 3/4 then ^ 



n-1 



-n+1 



\p{s)\ = 0{n 



(see Lemma 



I3.12P and X]rezP^(^) < oo so that, in view of (I3.38p . one has Var(i||Z}Z„|||,) = 0( 



n 



U H = 3/4 then ^ 
Lemma 13.121 p^(r 



n-1 
s=— n+1 
est 
|r| 



\pis)\ 



as r 



0(v/ri) and Zlr=-n+i = O(logn) (indeed, by 

oo) so that, in view of (13.39p . Var(|||DZ„||^) 



0(1/ log n). Finally, the desired conclusion follows from Theorem 13.61 



□ 



3.5 The method of (fourth) moments: exphcit estimates via in- 
terpolation 

It is clear that the combination of Theorem 13.61 and Theorem 13.71 provides a remarkable 
simplification of the method of moments and cumulants, as applied to the derivation of 
CLTs on a fixed Wiener chaos (further generalizations of these results, concerning in par- 
ticular multi-dimensional CLTs, are discussed in the forthcoming Section HI). In particular, 
one deduces from (I3.35p that, for a sequence of chaotic random variables with unit vari- 
ance, the speed of convergence to zero of the fourth cumulants E{Z^) — 3 also determines 
the speed of convergence in the Kolmogorov distance. 
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In this section, we shall state and prove a new upper bound, showing that, for a 
normalized chaotic sequence {Z„ : n ^ 1} converging in distribution to ~ -/^(0, 1), the 
convergence to zero of E{Z^) — E{N^) is always dominated by the speed of convergence of 
the square root of E{Z^) — E{N'^) = E{Z^) — 3. To do this, we shall apply a well-known 
Gaussian interpolation technique, which has been essentially introduced by Talagrand (see 
e.g. [HHj); note that a similar approach has recently been adopted in [22], in order to deduce 
a universal characterization of CLTs for sequences of homogeneous sums. 

Remark 3.13 1. In principle, one could deduce from the results of this section that, 
for every k ^ 3, the speed of convergence to zero of kth cumulant of Zn is always 
dominated by the speed of convergence of the fourth cumulant E{Z^) — 3. 

2. We recall that the explicit computation of moments and cumulants of chaotic random 
variables is often performed by means of a class of combinatorial devices, known as 
diagram formulae. This tools are not needed in our analysis, as we rather rely on 
multiplication formulae and integration by parts techniques from Malliavin calculus. 
See [30l Section 3] for a recent and self-contained introduction to moments, cumulants 
and diagram formulae. 

Proposition 3.14 Let q ^ 2 be an integer, and let Z be an element of the qth chaos Hq of 



X. Assume that Var(Z) = E{Z^) = 1, and let N 



\E{Z^) - E{N'')\ ^ Ck,q^/E{Z'^)-E{m] 
where the constant Ck,q is given by 



^(0, 1). Then, for all integer k ^ 3, 

(3.40) 



(A; - 1)2'=- 



3q 



i^^^ + (2.-5)^- 



2'=-2(A;-2)! 



Proof. Without loss of generality, we can assume that is independent of the underlying 
isonormal Gaussian process X. Fix an integer k ^ 3. By denoting \E'(t) = E [{y/l — tZ + 
ViN)''], t e [0,1], we have 



where the derivative \E'' is easily seen to exist on (0, 1), and moreover one has 



^ E[{VT^Z + ViN)''~^N] -^=E[{VT^.Z + ViNf~^Z]. 



By integrating by parts and by using the explicit expression of the Gaussian density, one 
infers that 

E[{VT^Z+VtNf'^N] = ElEliVT^z + ViNf'^N]^^^^ 

= {k-l)ViE\E[{Vr^z + ViN)''-^]^^^^ 
= {k-l)ViE[{VT^Z + ViN)''^^]. 
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Similarly, using this time (12.221) in order to perform the integration by parts and taking 
into account that {DZ, —DL^^Z)fj = because Z G Hg, we can write 



E[{VT^Z + ViN'f-^Z] = E 



E[{VT^tZ + Vtxf~^Z]^^^^ 

E[{Vr^Z + Vixf-^-\\DZ\\l] 



{k-l)Vl-tE 
{k - l)Vl -tE 



\x=N 



{VT^Z + ViN)''-'^^\\DZ\\l 



Hence, 

^ ' 2 
and consequently 



1 - -\\DZ\\% ] {VT^z + VtN) 



k-2 



\M;'(t)\^^^^-^jE 



\ 



E 



\DZ\\l 



By (13:271) and iKT2\\ . we have 

2" 



E 



\DZ\\l 



and E{N 



Using succesively {x + yY^ ^ ^ 2^^ ^{x"^^ ^ + ?/^^' ^Jx + y ^ v^+^, inequality (I2.12p 



{2k - 4)!/(2^-2(^ _ 2)!)^ can write 



^ 2^'"§(1 - t)t-i(2A; - 5)^"^ + 2^-5tf-i 



(2A;-4)! 
2'=-2(A;-2)! 



so that 



(2A;-5) 



{2k-A)\ 
2*^-2(A;-2)! 



Putting all these bounds together, one deduces the desired conclusion. 



□ 



4 Multidimensional case 

Here and for the rest of the section, we consider as given an isonormal Gaussian process 
{X{h) : h G over some real separable Hilbert space 9). 
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4.1 Main bounds 



We shall now present (without proof) a result taken from [23], concerning the Gaussian 
approximation of vectors of random variables that are differentiable in the Malliavin sense. 
We recall that the Wasserstein distance between the laws of two M'^- valued random vectors 
X and y, noted dy/{X,Y), is given by 

d^{X,Y) := sup \E[g{X)]-E[g{Y)]\, 

where indicates the class of all Lipschitz functions, that is, the collection of all functions 
g ■.R'^ such that 

II II \9{x) -9{y)\ ^ 
WgUip ■= sup — I II < oo 

(with II ■ ||]gd the usual Euclidian norm on M''). Also, we recall that the operator norm of a 
d X d matrix A over R is given by ||A||op := sup||^||^^^;^ ||/lx||igd. 

Note that, in the following statement, we require that the approximating Gaussian 
vector has a positive definite covariance matrix. 

Theorem 4.1 (See |23| ) Fix d ^ 2 and let C = {Cij)i^ij^d be a d x d positive definite 
matrix. Suppose that N ~ ^(0,C), and assume that Z = {Zi, . . . , Za) is a R'^-valued 
random vector such that E[Zi] = and Zi G D^'^ jor every i = 1, . . . , d. Then, 

d 

E[{a, - {DZ,, -DL-^Z,)^y]. 

In what follows, we shall use once again interpolation techniques in order to partially 
generalize Theorem 14.11 to the case where the approximating covariance matrix C is not 
necessarily positive definite. This additional difficulty forces us to work with functions that 
are smoother than the ones involved in the definition of the Wasserstein distance. To this 
end, we will adopt the following simplified notation: for every </) : M'^ — > M of class C^, we 
set 

dxidxj 

Theorem 4.2 (See |22| ) Fixd ^ 2, and let C = (Cjj)i^jj^rf be adxd covariance matrix. 
Suppose that N ~ ^/^(O, C) and that Z = {Zi, . . . , Za) is a R'^ -valued random vector such 
that E[Zi] = and Zi G D^'^ for every i = 1, . . . ,d. Then, for every (p ■.R'^ ^ R belonging 
to such that ||v^"||oo < oo, we have 

1 

\E[if{Z)] - E[ip{N)]\ ^ -y"\\^ Yl ^ O^M - {DZj, -DL-^Z,)^\] . (4.41) 



d^{Z,N) ^ II C 



I op 



\C 



1 1/2 
I op 



\ 



Woo = max sup 
ij=i,...,d^ggd 
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Proof. Without loss of generality, we assume that is independent of the underlying 
isonormal Gaussian process X. Let </) : M*^ ^ M be a C^-function such that ||v5"||oo < oo. 
For any t E [0, 1], set ^(t) = E[ip{Vl - tZ + v^iV)], so that 

\E[ip{Z)] - E[ip{N)]\ = 1^(1) - ^(0)| ^ / 

Jo 

We easily see that \i/ is differentiable on (0, 1) with 



2v^ 2y/T^ 



Z, 



By integrating by parts, we can write 

dip 



E 



dxi 



ElE 



{y/T^Z + ViN)Ni 



J \z=Z _ 



d 



J \z=Z _ 



By using (12.221) in order to perform the integration by parts, we can also write 



E 



^{VT^tz + ViN)Zi 



ElE 



dip 
dxi 



{VT^Z + Vix)Zi 



\x=N 



VT^tJ2E{E -^{VT^tZ + Vlx){DZ,,-DL~'Z,)^ 



\x=N ^ 



r Pfi 

^^E^ Q^\^^tZ + VtN){DZ,,-DL-'z:)^ 



Hence 



1 " (9^ 
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so that 



"'O i,j=l 

and the desired conclusion follows. 

We now aim at applying Theorem 14.21 to vectors of multiple stochastic integrals. 



□ 



Corollary 4.3 Fix integers d ^ 2 and 1 ^ gi ^ . . . ^ g^. Consider a vector Z = 
(Zi, . . . , Za) := (Jg,(/i), . . . , Jg,(/rf)) with U e ^1®^' for any 1 = 1..., d. Let N ~ C), 
with C = {Cij)i^ij^d a d X d covariance matrix. Then, for every : M"^ — >■ M belonging to 
such that ||¥'"||oo < oo, we have 



1 

\E[^{Z)]-E[v{N)]\^-y"Uj2^ 



Cij — —{DZj,DZi)f, 



(4.42) 



Proof. We have —L ^Zi = j- Zi so that the desired conclusion follows from (14.4ip . 



□ 



When one applies Corollary 14.31 in concrete situations, one can use the following result 
in order to evaluate the right-hand side of fl4.42p . 

Proposition 4.4 Let F = Ip{f) and G = Ig{g), with f G S)®^ and g e fl®^ (p,q ^ 1). 
Let a be a real constant. If p = q, one has the estimate: 



E 



a--{DF,DG), 



2 



r=l 



1 {2p - 2r)\{\\f 'S)p^r f\\ls>2r + \\g 0p^r g\\%Si2r) . 



p — 1 
r 



On the other hand, if p < q, one has that 

2 



E 



a-- {DF,DG)^ 



^ a + p\ ' 



p 



{q-p)\ 



-<>v\\g ®q-p g\\s^<i>2p 



2 



r=l 



p — 1 

r — 1 



A (P + g - 2r)!(||/ ®p-r f\\%mr + 11^ ®q-r ^|||«2.) . 



Remark 4.5 When bounding the right-hand side of (I4.42p . we see that it is sufficient to 
asses the quantities \\fi ®r fi\\f)'»2{qi-r) for alH = 1, . . . , d and r = 1, . . . , gj — 1 on the one 
hand, and E{ZiZj) = qi\{fi, fj)sj^ii for all i,j = 1, . . . ,d such that g^ = qj on the other 
hand. In particular, this fact allows to recover a result first proved by Peccati and Tudor 
in [31], namely that, for vectors of multiple stochastic integrals whose covariance matrix 
is converging, the componentwise convergence to a Gaussian distribution always implies 
joint convergence. 
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Proof of Proposition \4.4\ Without loss of generality, we can assume that ^ = 

where {A, s^') is a measurable space, and is a a-finite and non-atomic measure. Thus, 

we can write 

{DF.DG)^ = pq{Ip-i{f),I,-iig))^=pq [ Ip-i{f{-,t))l,.,{g{-,t))f^{dt) 

J A 



pq 



pAq—l 



r=0 



p — 1\ f q — 1 



p+q-2-2r 



{fi;t)^rgi;t))fi{dt) 



pAq-l , 

pq ^ r\{ 

r=0 ^ 



pAq 

pq^{r - 1)! 



p — 1\ f q — 1 



r 



Ip+q-2-2r{f®r+ig) 



r=l 



p-l\fq-l 
r — 1 1 \ r — 1 



Ip+q-2rif®rg)- 



It follows that 



E 



a - - (DF, DG) 



(4.43) 



\2 I ^2Y^P-1/ 



+ q-2r) 

2 



if p < q, 



{a-p\{f,g)^^r.r+p' ZViir " lV'{lZ\) i^P " 2r)!||/®,r7|||«,,,„,., if p = q. 
If r < p ^ q then 

||/®r5'll|®(p+9-2r) ^ \\f ®r g\\%<s(p+<i-2r) = if ®p-r f,g ®q-r g)sj<S2r 

^ 11/ ®p-r /II f,®2r II 5( (g)g„r fi'|U»2r 
1 



^ 2 ("•^ /llf,®2r + \\g ®q^rg\\ls>2r) . 



If r = p < g, then 



Wf^pg'' 



p y iif,®(9-p) 



^ 11/ ®p g 



P i/|lf,8)(9-p) 



If r = p = g, then /(Sips' = {f,g)f,e>p- By plugging these last expressions into (14.431) . we 
deduce immediately the desired conclusion. |-| 



4.2 Quadratic variation of fractional Brownian motion, continued 

In this section, we continue the example of Section 13. 4[ We still denote by i? a fractional 
Brownian motion with Hurst index HE (0,3/4]. We set 

Zn{t) = — J2 - ^k)' - 1], t^O, 
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where cr„ > is such that = 1. The following statement contains the multidi- 

mensional counterpart of Theorem 13.101 namely a bound associated with the convergence 
of the finite dimensional distributions of : t ^ 0} towards a standard Brownian 

motion. A similar result can be of course recovered from Theorem 14.11 - see again |23) . 

Theorem 4.6 Fix d ^ 1, and consider = to < ti < . . . < 1^. Let N ~ ^(0, Id). There 
exists a constant c (depending only on d, H and ti, . . . ,td) such that, for every n ^ 1; 



sup 



E 



l<i<d 



-E[^iN)] 



^ c X < 



9 IJ o 

n 2 



2/ if G (0,1] 



ifH 



where the supremum is taken over all C'^ -function ip : 



^ Vlogn 4 

such that ||v5"||oo ^ 1- 



Proof. We only make the proof for if < 3/4, the proof for ii = 3/4 being similar. Fix 
d ^ 1 and to = < ti < . . . < t^. In the sequel, c will denote a constant independent of 
n, which can differ from one line to another. First, see e.g. the proof of Theorem 13.101 
observe that 



with 



J n 



ii) 



Zn{ti) — Zniti^i^ 
^/ti ~ ti-l 

1 



\nU\-l 



in)-. 



i-1 



E IS 



k+iy 



In the proof of Theorem 13.101 it is shown that, for any fixed i G {l,...,d} and r G 
{l,...,g,-l}: 



1 



n 



2H- 



if ii G (0,i] 



(4.44) 



Moreover, when 1 ^ i < j ^ c/, we have, with p defined in (13.370 . 

[nU\-l \ntj\-l 

E E 



1 



a. 



[ntj\-[nti-i}-l 

Yl [(K-J-i-Oa(KJ 

|r| = [nij_ij-[ntij+l 



^ C 



[nti\ - [nti„ij - 1 



|r|>[nfj_ij-[nfij+l 



(K_iJ-r)V(K_iJ)]p2(r 

(4.45) 



as n — >• oo, 
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the last equality coming from fl3.38p and 

J2 p\r) = 0( ^ |r|^^-^) = 0(iV4^-3), as iV ^ oo. 

Finally, by combining (14.441) . (14.450 . Corollary 14.31 and Proposition 14.41 we obtain the 
desired conclusion. ^ 
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